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We examine the performance of a TPC as a y-ray telescope above the pair-creation threshold. The contributions to the photon 
angular resolution are studied and their dependence on energy is obtained. The effective area per detector unit mass for such a thin 
detector is the conversion mass attenuation coefficient. The differential sensitivity for the detection of a point-like source is then 
derived. Finally, the measurement of track momentum from deflections due to multiple scattering is optimized. 

These analytical results are exemplified numerically for a few sets of detector parameters. TPCs show an impressive improvement 
in sensitivity with respect to existing pair-creation-based telescopes in the [MeV - GeV] energy range, even with the modest detector 
parameters of this study. In addition, gas TPCs allow an improvement in angular resolution of about one order of magnitude. 
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1. Introduction 

High-energy astrophysics focuses on the study of com- 
plex objects such as active galactic nuclei (AGN), pulsars 
and gamma-ray bursts (GRB), the spectral energy distributions 
(SEDs) of which comprise contributions from a number of pro- 
cesses Q], including synchrotron emission, synchrotron self- 
Compton scattering (SSC), the inverse scattering of external 
thermal photons from the accretion disk (ESC), and tt° decays. 
These studies suffer from the existence of gaps in the sensitivity 
of the available instruments, such as that between the Compton- 
scattering-based telescopes, that are mostly sensitive below 1 
MeV, and the pair-based telescopes, that are mostly sensitive 
above 100 MeV. Pair-based telescopes also suffer from a poor 
angular resolution at low energy, making the analysis of dense 
regions of the sky difficult. 

In this paper we present a characterization of the use of a thin 
detector as an active target for cosmic y-ray detection in the 
pair-conversion regime, with the aim of improving the angular 
resolution and the sensitivity at low energy. The past missions 
and the presently operated Fermi/LAT and AGILE use thick de- 
tectors for which the photon conversion probability is close to 
unity and the effective area A e ff is the product of the geomet- 
rical area of the detector A and of the reconstruction efficiency 
€ r : A e ff = A x e r . In a thin detector for which the conversion 
probability is small, the effective area becomes proportional to 
the pair conversion mass attenuation coefficient H and to the 
detector mass M, A e gp = H x M x e r . In a thick detector, the 
optimization of the aspect ratio height/surface is therefore criti- 
cal, while for a thin detector it does not affect the effective area 
to first order. Note that in a thick detector, the various pos- 
sible conversion processes are in competition with each other, 
and Compton scattering prior to- or in place of- pair conversion 
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would become a problem at the lowest energy if low Z material 
were considered, an effect strongly reduced for a thin detector. 
A sample of previous works on the use of gas detectors for the 
detection of cosmic y rays above the pair-creation threshold, 
sometimes in relation with a Compton telescope, can be found 

inRefs. man. 

2. Test detector parameters 

We use a time projection chamber (TPC) as a test model of a 
thin detector. A TPC is a chunk of homogeneous matter, located 
inside a (say) uniform electric field E ([5 1, and a recent review 
[6]). Charged tracks crossing the TPC ionize atoms/molecules 
on their path, after which the ionization electrons drift to the 
anode plane where (in a gas TPC) they are amplified, and col- 
lected. The collecting plane can be segmented so as to provide 
two coordinates x, y of the location of the ionization. The third 
coordinate z is provided by the drift time t, z - vj x t. The 
drift velocity vj ranges from (9(mm/|j.s) at saturation for liquids 
and solids, to Vd = <9(cm/|j.s) for gases when used with an ap- 
propriate quencher, a multi-atomic molecule on which drifting 
electrons collide and are cooled down. The size of the electron 
"cloud" increases during the drift due to diffusion, which places 
a limit to the useable drift length. The typical value of the diffu- 
sion coefficient is (9(100 - 200 \im/ Vcm) for the three phases 
considered here. 

Note that amplification of the signal is easily performed in 
a gas, but not in a liquid. The smallness of the collected sig- 
nal * 4000 e~ /mm in liquid argon (lAr) and the limitation of 
the electrical power available in space for the digitization elec- 
tronics would make precision tracking difficult. To alleviate 
this limitation, double-phase systems have been developed: the 
electrons are extracted from a liquid TPC into a gas in which 
they are amplified efficiently [7 8 1. But for a use in space, the 
instability of the gas/liquid interface would be a problem. A 
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gas/solid system might be considered such as solid argon and a 
neon-based gas. The density and electron parameters of noble- 
gas solid and liquid are close enough that we neglect the dif- 
ference and we use the liquid as representing the two dense 
phase^j 

These analytical results are exemplified numerically for a few 
sets of detector parameters : three noble gases: neon, argon and 
xenon; three "densities": gas at 1 and 10 bar, and liquid. We 
chose on purpose rather modest detector parameters, a detector 
mass M — 10 kg and a tracking length of L — 30 cm for gas, 
M — 100 kg and L = 10 cm for liquid, except in the case where 
the radiation length Xq would be shorter than L: in that case 
the propagation of an electron would be affected by radiation 
after a path length of the order of Xq, and we use L = Xq. The 
tracking length L is obviously related to the detector thickness, 
but they can be different, as in the case of liquid xenon. The 
sensitive mass for a gas TPC is chosen to be smaller due to the 
volume limitation of space missions. 

We use typical values for the longitudinal sampling of the 
TPC, / = 1mm, and of the point space resolution, cr = 0.1mm. 

3. Angular resolution 

We consider a photon with energy E converting to an 
electron-positron pair in the field of a charged particle of the 
detector. The conversion is said to be "nuclear" in the case of 
an ion, yZ — > e + <?~Z, and "triplet" in the case of an electron 
ye~ — > e + e~e~. From momentum conservation: 
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we can identify several contributions to the photon angular res- 
olution: (1) the single track angular resolution; (2) the fact that, 
in the case of nuclear conversion, the recoil momentum q of 
the ion, of the order of 1 MeV/c produces a path length that is 
too small to allow a measurement of q; (3) the resolution of the 
norm of the momentum, named here the momentum resolution, 
for each of the hard tracks. We examine these three contribu- 
tions in the following. 

3.1. Single track angular resolution 

The basis of the understanding of tracking in the presence of 
multiple scattering was settled by Gluckstern flOl . We make 
use of the similar results derived by Innes for optimal fits [ 11 1. 
We use an approximation of the multiple scattering angle 9q un- 
dergone by a particle of momentum p crossing a slab of matter 
with thickness x lfl2"ll : 



scattering can be neglected and the detector resolution domi- 
nates 1111: 
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where N is the number of samplings, N - L/l. In the expres- 
sions of erg, the subscript t refers to a track and H and L to the 
high and low track momentum, respectively. o- etH is indepen- 
dent of the track momentum and improves like Zr 3 ^ 2 at given 
sampling I. On the contrary, multiple scattering dominates at 
low momentum; from the results of Ref. flTT l and using the 
approximation of eq. we obtain: 
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The expression does not depend on the tracking length. This 
is because the end of the track contributes no information about 
the particle's direction at the detector entrance once there has 
been sufficient scattering IfTTI . The p~ 3 ^ dependence on mo- 
mentum obtained here is comparable to the £ -0 - 78 dependence 
on energy parametrized by the Fermi Collaboration [13]. The 
momentum limit between the two regimes is at p — pu m , with: 
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where a p = (2/8 4 x 9) 
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0.038. 



3.1.1. From tracks to photon 

Combining the measurements of the directions of the elec- 
tron and of the positron in the small-angle approximation, the 
direction of the reconstructed photon (with respect to the (un- 
known) true incoming direction) is 6L = r9 x+ + (1 - r)6 x _, where 
r is the fraction of the energy that is carried away by the positron 
and 6 X+ and 6 X - are the angles of the positron and of the elec- 
tron, respectively. We then compute the photon angular resolu- 
tion from the track angular resolution, eqs. ([3j,Q. 
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where po = 13.6MeV/c and we have neglected the negative 
logarithmic correction factor. At high momentum, multiple 



191. 



Note that liquid neon does not allow electron to drift, but solid neon does 



Figure 1 : Track-to-photon angular resolution factors as a function of the frac- 
tion r of the photon energy carried away by the positron. 



In the multiple scattering dominated regime we obtain cr gy — 
CTet J yjr + Vl - r and in the high energy regime: cr gy = 

o~et yr 2 + (1 - r} 2 , in the expressions of which the track mo- 
mentum p is replaced by the photon energy E. The photon 
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angular resolution is therefore obtained from the single track 
angular resolution by applying a correction factor that is close 
to unity (Fig[T]l and that we neglect in the following. 
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Figure 2: (Color online). Dependence of the RMS photon angular resolution 
on photon energy, in argon for various densities (left) and in 10 b gas for various 
gases (right). 



Assuming a quadratic sum of <x e ,# and of cr 9fL in the inter- 
mediate energy range, we can examine the energy dependence 
of the photon angular resolution due to single track angular res- 
olution (Fig. El. Less dense TPCs have a better resolution. 
The difference in the high-energy asymptotic values, between 
gas and liquid, is due to the shorter tracking length used here, 
10 cm vs 30 cm. 



that way. The estimations of q^ from the ID q distribution 1 14| 
and from the 5D Bethe-Heitler differential cross section lfT31 
are found to be compatible with each other which constitutes a 
cross check (Fig. ffl. Also of interest is the momentum at half- 
maximum q\j2, on the decreasing slope above the maximum. 
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Figure 4: (Color online). Left: 68 % "containment" value q^g of the recoil 
momentum q, as a function of incoming photon energy E. 'J' stands for the q 
distribution by Jost et al. (squares) 1141 . while 'BH' stands for the 5D Bethe- 
Heitler differential cross-section (bullets) 1151 . Values qM of q at maximum 
(up triangles) and at half maximum qi/2 (down triangles). The solid symbols 
corresponds to the absence of screening, and the open symbols to screening 
parametrized by the Mott form factor [ IT]. Right: 68 % "containment" space 
angle induced by the non-observation of the recoiling nucleus. The upper line 
corresponds to 8 = 1.5 rad[£/l MeV] -5 ' 4 , the lower lines to qi/2/E and qu/E. 



3.1.2. Recoiling ion 

In the case of nuclear conversion, the recoil momentum goes 
undetected, unless a very low pressure TPC is used. The recoil 
is almost transverse to the photon direction, and the contribu- 
tion to the photon angular resolution is therefore as q/E. The q 
distribution for nuclear conversion has been obtained by Jost et 
al. lfT4l Q by the integration of the 5D Bethe-Heitler differen- 
tial cross section fl5l . It has a very wide spectrum that extends 
from 2m 2 /E to E at high energy, where m is the electron mass. 
The spectrum peaks at a value of the recoil momentum q^ that 
decreases with E (Fig. 15). 
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Figure 3: (Color online). Distribution of the recoil momentum for several in- 
coming photon energies, from Ref. 1 14|, in the approximation of no screening, 
normalized to unity at their maximum. 



The 68 % containment value q^ of q is of particular interest 
as the point spread function of y telescopes is often expressed in 



Under the approximation of no screening of the field of the 
nucleus by the surrounding electrons, we find qM and q\/2 to 
vary like w 2.24MeV/(c£') and « 9.28 MeV/(c£), respec- 
tively, while ^68 has a much milder dependence on E, as q^ » 
(1.5 MeV /c)E~ l l 4 . Therefore, at high energy, a low-g core de- 
velops that is narrower with increasing E, but that contains a 
fraction of events which decreases with E. 







No screening 


Screening 






(low energy) 


(high energy) 




q e 


q e 


<?68 


68 % containment 


oc £~ 1/4 oc £- 5/4 


constant oc 1/E 




at maximum 


oc 1/E oc 1/E 2 


constant oc 1/E 




at half max 


oc 1/E oc 1/E 2 


constant oc 1/E 



2 With Borsellino's correction 1 16| applied. 



Table 1: Dependence on photon energy of q(,g, qM, and 171/2 and related angular 
shifts from 8 a q/E. 



The dependence of qM, qi/2 and q^ on photon energy are 
represented in Fig. |4] and tabulated in Table[T] Closed symbols 
refer to the absence of screening. The screening of the field 
of the nucleus by the surrounding electrons can be described 
by multiplying the q distribution [14] by (1 - F(q)) 2 , where 
F(q) is the Mott atomic form factor F(q) - l+ ^/ qo y_ , with qo = 
111Z~'/ 3 x m [17]. When screening is taken into account, the 
68% containment value saturates at high energy to a value that 
has a mild Z dependence, and that is close to 0.4 MeV/c, while 
the distribution maximum q M saturates to a much smaller value 
close to 2 x qo. The induced photon angle shift is obtained 
simply by 9 = q/E (Fig. |4]right). 



3.2. Contribution of track momentum resolution 

The contribution of the track momentum resolution Ap in the 
reconstruction of the photon momentum is estimated using a 
generator of the 5D Bethe-Heitler [ 15 1 differential cross section, 
to which a momentum spread assumed to be of 10 % is added. 
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Figure 5: (Color online). 68% containment angle shift due to a track momen- 
tum resolution of Ap/p = 10%, from a generation of the 5D Bethe-Heitler 1 15 1 
differential cross section. 



The dependence is found to be 6 m (Aplp)Q.% md/E (Fig. |5j. 
This l/E dependence was expected given the l/E dependence 
of the opening angle 1 18 1. 

3.3. Angular resolution: a Summary 

The results obtained so far are summarized in Fig. [6](we have 
neglected the small difference between the RMS and the 68% 
containment values). 



A ten-fold improvement of the angular resolution with re- 
spect to the performance of the Fermi/LAT (13 ] is within reach 
with a gas TPC. A liquid xenon TPC shows a poor performance 
here, in contrast with its use as a Compton telescope| 19]. 



4. Effective area 

Assuming 100% reconstruction efficiency, which is reason- 
able for a TPC, we obtain the dependence on energy of A e g- 
(Fig. [7]i from the tabulation of the mass attenuation coefficient 
inRef. - 
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Figure 7: (Color online). Dependence on energy of the effective area per unit 
mass assuming 100% reconstruction efficiency 1201 . 
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Figure 6: (Color online). Various contributions to the photon angular resolu- 
tion. 



In the case of triplet conversion, the effective area is almost 
independent of Z which is not surprising, as the number of elec- 
trons per unit mass is almost a constant of Nature. 



5. Differential sensitivity 

From the above we compute the sensitivity to the detection 
of a faint high-latitude point-like source. We mimic the compu- 
tation of the Fermi/LAT sensitivity 1 13 1 with a simpler method. 

We use four bins per energy decade (AE = lO^ 4 ^ « 
0.78E), a n = 5 standard deviation significance, a T = 3 year 
observation time, an rj = 17% exposure fraction, and a mini- 
mum signal event number S > 10. The sensitivity is computed 
inside a <x e 68% -containment angle, and "against" the extra- 
galactic y-ray background ED : we therefore compare our re- 
sults to Fermi's differential sensitivity at 90° galactic latitude 

m. 



• The single track angular resolution dominates; 

• The missing momentum recoil contributes significantly at 
low energy, especially for low density TPC; 

• The effect of a single track momentum resolution at a level 
of 10 % as assumed here is negligible. 



5.7. Gaussian statistics 

Using Gaussian statistics, n is simply obtained from the num- 
bers S and B of signal and background events, respectively, 
from S = n V#. We compute these numbers of events as : 



S = T x i] x A eff (£) x /()(£) X£XA£ 



(6) 
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B = Txj]X A eff (£) x no-g F B {E)dE 



(7) 



where Fb is the background flux, /rj the signal intensity, and 
e = 0.68, the efficiency of the angle cut. The sensitivity s ex- 
pressed as the minimum detectable signal intensity, multiplied 
by E 2 , that is s — E 2 x /q, is therefore: 
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• In the spatial resolution dominated regime (high momen- 

. % °> 1 / P f u PoA 3/2 , . . 
turn): — = — , where p m - , and A is 

P y/2N\Pm) o-a/Xo 

the length over which the deflection is measured, which 

must be larger than the sampling of the TPC, A > /. 

Approximating the full expression by the sum of the two 
asymptotes, we obtain: 
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and A e fj oc H 
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that is, asymptotically, A eff oc 1 /X . These dependences cancel the minimum of which is obtained for: 
partially in the expression of s oc erg/ jAlfr, so that s oc X. 



1/8. 
_0 • 

the Z dependence of the sensitivity as described by eq. <JSJ is 
extremely small. 

5.2. Poisson statistics 

Using Poisson statistics, and following the method used by 
Ref. [22] (A recent derivation of eq. |9]l can be found in Ref. 
Il23ll eq. (96)), the significance n becomes such that: 



n 2 = 2B[(l +5/B)ln(l +S/B)-S/B], 



(9) 



which reduces to n 2 = S 2 /B in the background-dominated 
regime S/B <s: 1, as expected. We solve eq. (|9]l for S, with 
n = 5. 

Figure [8] shows the dependence on energy of the differen- 
tial sensitivity for the nine test cases of this study compared 
to that of the Fermi/LAT iTPJl and of the Compton telescope 
COMPTEL EH. On the left plot, an ideal case for which the 
missing recoil momentum could be neglected in the expression 
of the angular resolution, we see that a gas TPC would clearly 
have a better performance than a liquid TPC: the better angular 
resolution would outweigh the smaller sensitive mass (10 kg 
compared to 100 kg). In reality, the effect of the missing re- 
coil momentum on the sensitivity is visible at low energy (right 
plot). 

Together with Compton telescope projects, which aim at sen- 
sitivities of the order of 10~ 5 MeV/(cm 2 s) close to 1 MeV 1251 . 
there is good hope to fill the sensitivity gap between the energy 
ranges of the Compton and pair telescopes, and even to envis- 
age enough overlap for a cross-calibration. 

6. Track momentum from multiple measurement of multi- 
ple scattering 

In this section we examine the potential of a TPC to deter- 
mine the momentum of the conversion electrons from the de- 
flection due to multiple scattering: since the average deflec- 
tion angle #o is proportional to 1/p, each deflection provides 
a momentum measurement l26]| . Bolton describes the relative 
momentum resolution cr p /p as asymptotes of a "universal func- 
tion" E71: 

• In the multiple scattering dominated regime (low momen- 

CT p 1 

turn): — = , where N is the number of samplings. 
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As the momentum is not known a priori, some iteration will 
be needed. The value of <r p /p for that optimal set is: 
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and C s 5 1/6 + 5~ 5/6 ~ 1.57. The method is usable at low 
energy, below 100 MeV (Fig. [9}) 
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Figure 9: (Color online). Relative track momentum resolution for optimal 
sampling A as a function of track momentum. (The low momentum asymptote 
for the liquids is due to the A > / limit). 



7. Summary 

We examined the performance of a TPC as a y-ray telescope 
above the pair-creation threshold. Analytical expressions are 
obtained for the various contributions to the photon angular 
resolution and the track momentum resolution. These analyt- 
ical results are exemplified numerically. Even with the modest 
value of the parameters used in the present study, a L — 30 cm 
tracking length and M — 10 kg sensitive mass for gases and 




E (MeV) E (MeV) 

Figure 8: (Color online). Dependence of the differential sensitivity on photon energy compared to the 90° galactic latitude performance of the Fermi/LAT [ 13] and 
of the Compton telescope COMPTEL [24]. Left : the ideal case in which the missing recoil momentum would be neglected in the expression of erg. Right : the real 
case, with the full expression of erg. 



L = 10 cm tracking length and M = 100 kg for liquids, TPCs 
show huge potential in the MeV-GeV range. 

The angular resolution is found to have an E~ 3 ^ 4 energy de- 
pendence in the multiple scattering dominated energy range, 
except for gas TPC at low energy, where the missing recoil mo- 
mentum dominates with a steeper dependence close to £~ 5 ^ 4 . 
An improvement of a factor better than ten in angular reso- 
lution with respect to the present tungsten-slab-converter tele- 
scopes is within reach for a gas TPC. This provides at least 
two orders of magnitude in background rejection, not to men- 
tion the easy albedo rejection and the immunity to cosmic-ray 
pile-up thanks to the detailed imaging provided by a TPC. The 
improved angular resolution translates also into a dramatic im- 
provement of the differential sensitivity that will bridge up with 
a future Compton mission at a couple of MeV at a level of about 
10- 5 MeV/(cm 2 s). 

Track momentum can be measured from the deflections due 
to multiple scattering across the detector. After optimization of 
the sampling size, the momentum relative precision is found to 
vary like p lli . The method can be used for momenta smaller 
than 100 MeV /c. Above that limit, the photon energy measure- 
ment must be performed by an additional device. A TPC might 
well be a first layer of an otherwise high-energy (0. 1 Ge V - 3 
TeV) range project such as Gamma-400 |28 1. 

An R&D program is in progress to characterize the use of a 
TPC as a y-ray telescope and as a y-ray polarimeter ||29l[30l . 

A number of approximations have been made in this work so 
as to obtain scaling laws in the form of analytical expressions 
that are easily handled. To name a few, we have neglected the 
logarithmic correction of multiple scattering, the loss of energy 



of the tracks inside the detector, their emission of radiation that 
disturbs tracking at high momentum, especially for high-Z de- 
tectors and/or dense TPC - such as liquid xenon, and trigger 
inefficiency. These limitations would be addressed using a full 
simulation, which is beyond the scope of the present study. 
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